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1 Introduction 



The quantum version of the inverse scattering method paved the way for the discovery and 
solution of new exactly solvable two-dimensional lattice models. Of particular interest are 
non-homogeneous vertex models whose transfer matrix are constructed by mixing different 
local vertex operators. These operators, known as Lax /^-operators, define the local structure 
of Boltzmann weights of the system. A sufficient condition for integrability is that all the 
/^-operators should satisfy the Yang-Baxter equation with the same invertible i?-matrix. More 
precisely we have |IJ, 

R(X - n)£ Ai (\) <g> C M (fj) = C M (ji) ® C Ai (X)R(X - //) (1) 

The auxiliary space A corresponds to the horizontal degrees of freedom of the vertex model 
in the square lattice. The operator C A i(X) is a matrix in the auxiliary space A and its matrix 

L 

elements are operators on the quantum space <8>Vi, where Vi represents the vertical space of 

i=l 

states and i the sites of a one-dimensional lattice of size L. The tensor product in formula (1) is 
taken with respect the auxiliary space and A is a spectral parameter. One way of producing a 
mixed vertex model is by choosing ^-operators intertwining between different representations 
Vi of a given underlying algebra. This approach has first been used by Andrei and Johannesson 
PI to study the Heisenberg model in the presence of an impurity of spin-S and by De Vega 
and Woynarovich || to construct alternating Heisenberg spin chains. Subsequently, several 
papers have discussed physical properties of the latter models M, [5], as well as considered 
generalizations to include other Lie algebras [0, §] and superalgebras || [10], fTT^ p^| . 



In this paper, which we are pleased to dedicate to James McGuire, we construct and solve 
mixed vertex models whose /^-operators can be expressed in terms of the generators of the 



braid-monoid algebra |]I3 |. In particular, we show that the recent results by Abad and Rios 
|| [H| and by Links and Foerster for the SU(3) and gl(2\l) algebras can be reobtained 
from such algebraic approach. The novel feature as compared to these works is that we are 
able to diagonalize the corresponding transfer matrix by using a standard variant of the nested 
Bethe ansatz approach. Furthermore, this algebraic approach allows us to derive extensions 
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to the SU(N) and Sl(N\M) algebras in a more direct way. 

This paper is organized as follows. In section 2 we recall the basics of the braid-monoid 
algebra and show how it produces two different /^-operators satisfying the Yang-Baxter equa- 
tion (1). In section 3 a mixed vertex model based on the SU(N) algebra is diagonalized by 
the algebraic Bethe ansatz method. 



2 Braid- monoid L-operators 

It is well known that the braid algebra produces the simplest rational i?-matrix solution of the 
Yang-Baxter equation. In this case the braid operator becomes a generator of the symmetric 
group and the i?-matrix is given by 

R(X) — Ii + Xh (2) 

where Jj is the identity and h\ is the braid operator acting on the sites % and i + 1 of a one- 
dimensional chain. Here we choose the braid operator as the graded permutation between N 
bosonic and M fermionic degrees of freedom [p~L| 



N+M 

b*= E (-l) P(aM6) 4 ® el? (3) 

a,b=l 

where p(a) is the Grassmann parity of the a-th degree of freedom, assuming values p(a) = 
for bosons and p(a) = 1 for fermions. 

The i?-matrix (2,3) has a null Grassmann parity, and consequently can produce a vertex 
operator Cm{X) satisfying either the Yang-Baxter equation or its graded version |L4j]. In the 
latter case, the tensor product in formula (1) should be taken in the graded sense(supertensor 
product) |TJ|. In the latter case the associated ^-operator is 



C M {\) = XI t + b l (4) 

The next step is to search for extra /^-operators which should satisfy equation (1) with the 
i?-matrix (2,3). As we shall see below, this is possible when we enlarge the braid algebra by 
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including a Temperley-Lieb operator Ei. This operator satisfies the following relations 

Ef = qEi, EiE i±1 Ei = E h E { Ej = EjEt \i - j\ > 2 (5) 

where q is a c-number. It turns out that the braid operator (3) together with the monoid Ei 



close the braid-monoid algebra [|13j at its degenerated point {bf = Jj). The extra relations 



between bi and Ei closing the degenerated braid-monoid algebra are (see e.g flT5|) 

hiEi = Ejbj = tEi 

(6) 

Eibi ±1 bi = b i ± 1 b i E i ± 1 = EiE i±1 

where the constant t assumes the values ±1. 

Now we have to solve the Yang-Baxter equation (1) with the i?-matrix (2,3) assuming the 
following general ansatz for the ^-operator 

A/u(A) = /(A)ij + g(\)bi + h{X)Ei (7) 

where /(A), g(X) and h(X) are functions to be determined as follows. Substituting this ansatz 
in equation (1), and taking into account the braid-monoid relations, we find two classes of 
solutions. The first one has h(X) = and clearly corresponds to the standard solution already 
given in equation (4). The second one is new and is giving by g(X) = 0. We find that the new 
/^-operator, after normalizing the solution by function h(X), is given by 

C Ai {X)=t{X + r ] )I i -E i (8) 

where rj is an arbitrary constant. This constant can be fixed imposing unitary property, i.e 
C^i(X)£^i(— X) ~ I{. By using this property and the first equation (5) we find 

" = l < 9 > 

After having found two distincts /^-operators which satisfy the Yang-Baxter algebra with 
the same i?-matrix, the construction of an integrable mixed vertex model becomes standard 
H §| . The monodromy matrix of a vertex model mixing L\ operators of type £/y(A) and L2 
operators of type Cm{X) is written as 

T L ^(X) = C ALi+L2 (X)£al 1+ l 2 ~i(X) ■ ■ -C M (X) (10) 



where £^«(A) is defined by 

\ C Ai (X) ifie {&,•••, m , 

I Adi(A) otherwise 

and the partition {/3i, • • • , (3^} denotes a set of integer indices assuming values in the interval 
1 < oci < Li + L 2 . Although the integrability does not depend on how we choose such partition, 
the construction of local conserved charges commuting with the respective transfer matrix does. 
One interesting case is when the number of operators £^i(A) and Cj^{X) are equally distributed 
(Li = L2 = L) in an alternating way in the monodromy matrix ||. In this case, the first non- 
trivial charge, known as Hamiltonian, is given in terms of nearest neighbor and next-to-nearest 
neighbor interactions. More specifically, the expression for the Hamiltonian in the absence of 
fermionic degrees of freedom (M = 0) is 

2L 2L 
H ^ ] £m— l,m 

(0)+ 4-2,n-l (0)£ n ,n-l (0)P„_2,„ (12) 

mGodd iGeven 

~ ~ 2 

where in the computations it was essential to use the unitary property £^i(0)£^(0) = ^-/j 
at the regular point A = 0. We also recall that Pjj denotes permutation between sites i and 
j (equation (3) with M — 0). 

We close this section discussing explicit representations for the monoid E%. Such represen- 
tations can be found in terms of the invariants of the superalgebra Osp(N\2M), where N and 
2M are the number of bosonic and fermionic degrees of freedom, respectively. Here we shall 
consider a representation that respects the U(l) invariance, which will be very useful in Bethe 
ansatz analysis. Following ref. |15| the monoid is written as 

N+2M 

Ei= J2 o^a^e^ ® ejj 1 (13) 

a,b,c,d=l 

and the matrix a has the following block anti-diagonal structure 



a 



^ OnxM OnxM Inxn ^ 

OmxM xM OmxN 

\ —Imxm OmxM OmxN J 



(14) 



where Xfc lX fc 2 and Ok lX k 2 are the anti-diagonal and the null k\ x k 2 matrices, respectively. We 
also recall that for t = 1 the sequence of grading is /i • • • fubi ■ ■ ■ &at/m+i • • • /2M and the 
Temperley-Lieb parameter q is the difference between the number of bosonic and fermionic 
degrees of freedom 

q = N-2M (15) 

Finally, we remark that new £41 (A) operators are obtained only when N + 2M > 3. 
Indeed, for the special cases N = 2, M = and iV = 0, M = 1 it is possible to verify that 
such operator has the structure of the 6- vertex model, which is precisely the same of C^(X), 
modulo trivial phases and scaling. In the cases N = 3,M = and N = 1,M = 1 we reproduce, 
after a canonical transformation, the ^-operators used recently in the literature || [10], [LT| to 
construct mixed 577(3) and t— J models. As we shall see in next section, however, an important 
advantage of our approach is that representation (14) is the appropriate one to allows us to 
perform standard nested Bethe Ansatz diagonalization. 

3 Bethe ansatz diagonalization 

In this section we look at the problem of diagonalization of the transfer matrix T Ll,Z/2 (A) = 
Tr A [T Ll ' L2 (X)], namely 

T Ll ' i2 (A) |$) = A(A) |$) (16) 

by means of the quantum inverse scattering method. 

For sake of simplicity we restrict ourselves to the case of mixed vertex models in the absence 
of fermionic degrees of freedom. In this case the ^-operators C Ai (A) and C Ai (A) are given by 
formulae (4) and (8) with M = 0, t = 1 and r\ = N/2. An important object in this framework 
is the reference state |0) we should start with in order to construct the full Hilbert space |<&). 
The structure of the ^-operators suggests us to take the standard ferromagnetic pseudovacuum 
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as our reference state, i.e 



Li+l 2 

10)= n ®io>*. i«>i 
i=i 





w 



(17) 



N 



where the index N represents the length of the vectors |0)^ It turns out that this state is an 
exact eigenvector of the transfer matrix, since both operators £4$ (A) and jCj^(X) satisfy the 
following important triangular properties 

fa(X)\0) i * * ... * \ 
6(A) 10), ... 



and 



V 

/6(A) 













6(A) |0). 





b W\o)J NxN 
* \ 



(1* 



(19) 



S(A)|0)J 



NxN 



V . 

where the symbol * stands for non-null values that are not necessary to evaluate in this algebraic 
approach. The functions a(A), 6(A), a(A) and 6(A) are obtained directly from expressions (8,9), 
and they are given by 



a(A) = A + l, 6(A) = A, a(A) = A + y-l, 6(A) = A + j 



(20) 



To make further progress we have to write an appropriate ansatz for the monodromy matrix 
T Ll,L:2 (A) in the auxiliary space A. The triangular properties (18,19) suggest us to seek for 



T 



(A) 



(21) 



standard structure used in nested Bethe ansatz diagonalization of SU(N) vertex models 

A(X) Bi(X) 
Q(A) A,(A) /ArxJV 

where i, j — 1, - • • , N—l. As a consequence of properties (18,19) we derive how the monodromy 
matrix elements act on the reference state. The fields -Bj(A) play the role of creation operators 



while Cj(A) are annihilation fields, i.e Cj(A) |0) = 0. Furthermore, the action of the "diagonal" 
operators A(X) and -Djj(A) are given by 

A(X)\0) = [a(X)}^[b(X)}^\0) (22) 

D U (X) |0) = [6(A)] Ll [&(A)] L2 |0) for i ^ N - 1, ZW^A) |0) = [6(A)] Ll [S(A)] ia |0> (23) 

A,i(A) |0) = for i ^ j and j=£ N- 1, A,jv-i(A) |0) ^ for i ^ JV - 1 (24) 

We observe that although matrix Dij(X) |0) is non-diagonal it is up triangular, which is 
an important property to carry on higher level Bethe ansatz analysis. In order to construct 
other eigenvectors we need to use the commutation relations between the monodromy matrix 
elements which are obtained by extending the Yang-Baxter relation to the monodromy matrix 
ansatz (21). Due to the structure of the i?-matrix, the commutation rules are the same of that 



already known for isotropic SU(N) models |1~6|, and the most useful relations for subsequent 
derivations are 

A(X)B^) = ^^B^)A{X) - ^±-^B t (X)A{fi) (25) 

D ij (X)B k (fM) = -^L_BMD iq {X)r^\X - n)f q - —L—B^DM (26) 

B % {X)B 3 {^) = B p ^)B q {X)r^\X - p)% (27) 

where ^^•'(A)^ are the elements of the i?-matrix Jj + Xbi on the subspace (N — 1) x [N — 1). 
The eigenvectors are given in terms of the following linear combination 



(Ai 1) 7 • • • , A«)) = B ai (X?) • • • B ami (A«)^-- ai (28) 

where the components jF am i " ai are going to be determined a posteriori. 

By earring on the diagonal fields A(X) and Du(X) over the above mi-particle state we 
generate the so-called wanted and unwanted terms. The wanted terms are those proportional 
to $ mi (Ai 1 ' ) , • • • , A^)) and they contribute directly to the eigenvalue A Ll ' L2 (A, {Aj- }). These 
terms are easily obtained by keeping only the first term of the commutation rules (25,26) 
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each time we turn A(X) and Da(X) over one of the B ai (X^) component. The result of this 
computations leads us to the following expression 

mi a(Aj 1) - A) 
S &(A? } - A) 



r tl - £a (A)|$ mi (Af ) > ...,A«)) 



•(A)]* [KA)] L2 II 777(1) — 77 |<*V (A?\ • ■ • , AS)) 
i=i 6(A) - A) 

mi 

+n 



,=i 6(A - A? } ) 
+unwanted terms 



where T^(A, {Xj }) is the transfer matrix of the following inhomogeneous auxiliary vertex 
model 



T (1) (A, {aS 1 )})^.:- = r«(A - A^/ 1 ^ - A«)^ ■ ■ -r«(A - A«)£^IU* (A) |0) 

(30) 

The unwanted terms arise when one of the variables A^ of the mi-particle state is exchanged 
with the spectral parameter A. It is known [16| how to collect these in a close form, thanks to 
the commutation rule (27) which makes possible to relate different ordered multiparticle states. 
We find that the unwanted terms of kind B ai (X^) ■ ■ ■ B a .(X) ■ ■ ■ B ami (A££) are cancelled out 
provided we impose further restriction to the mi-particle state rapidities A| , namely 

a(X^-X^) 



a(X 



6(A 



n & (a 

r( 1 )(A = A| 1 \{A«})^^.^ li 



A i ^7777i) 



6(A ■ J - A 



(i) 



_jra.m 1 —ai 



m 1 



(31) 



Now it becomes necessary to introduce a second Bethe ansatz in order to diagonalize the 
auxiliary transfer matrix T^(X, {A^ }). The only difference as compare to standard cases 
|jnS|| is the presence of the "gauge" matrix g ab = D a b(X) |0). It turns out that this problem is 
still integrable since the tensor product g <S> g commutes with the auxiliary i?-matrix r^'(X) 
[](see e.g ref. P7| ). From equations (23,24) we also note that this gauge does not spoil the 
triangular form of the monodromy matrix associated to T"W(A, {A| }) when it acts on the 

1 This occurs because the off-diagonal elements D^jv-^A) belongs to a commutative ring. 



S 



usual ferromagnetic state, 



i=i 







(32) 



N-l 



By defining A^^A, {A J - 1 ' ) }) as the eigenvalue of the auxiliary transfer matrix T"W(A, {A^}) 



i.e 



T( 1 )(A,{Al 1) })a 1 i :. fe a ': i 1 ^ mi --- bl = A (1) (A,{Af ) })^- ai 
we derive from equation (29) that the eigenvalue of T Ll ' L2 (A) is given by 



(33) 



A^(A,{A«}) = [a(A)]^[KA)] La n ? ,m , +11 



1 



jy-AWCA^AW}) (34) 



i=1 6(A 4 (1) -A) ti 6(A - A? } ) 
and the nested Bethe ansatz equations (31) become 



a(A 



(i)> 



6(A 



(i)> 



9 K^-^ffp, -A«>(A-A{'>,{^}), i 



6(Af } - A 



(35) 



In order to find the auxiliary eigenvalue A^(A = A^ , {A^}) we have to introduce a 



new set of variables {Ai , • • ■ , A^H which parametrize the eigenvectors of TW(A, {A| }). The 
structure of the commutations rules as well as the eigenvector ansatz (28) remains basically 
the same, and the expression for A«(A = Af^A^}) will again depend on another auxiliary 
inhomogeneous vertex model having (N — 2) states per link. We repeat this procedure until 
we reach the (N — 2)th step, where the auxiliary problem becomes of 6-vertex type. Since this 
nesting approach is well known in the literature [jlB] , we here only present our final results. 
The eigenvalue of the transfer matrix T Ll ' L2 (A) is given by 



A L " fa (A;{A< 1 >},...,{A<'"- 1 >}) 



+P>(A)] Ll |?>(A)] 12 H 1177^7(7)7 It ' 



K "i b(X - A*") b(Xj 

mjv " 1 n(\ — \^ N ~^\ 

j=i 6(A-A ') 



(i+i) 



A) 
(36) 



while the nested Bethe ansatz equations are given by 



(1)> 



(i)a 



n a(A« - Af) ™ a(Af - AM) 
a(A« - AS 1 )) M KAf - A?>) 



(37) 



a(A 



(0 



A<«) 



3=1 



6(A 



Af x) ) £i a(Aj 



P+i) 



Af } ) 



(3? 



a(A 



(iV-l)x 



6(A 



L 2 



mjv-i 

n 



a(A 



(jV-1) 



A 



(iV-l)N mjv _ 2 



a(A 



(iV-l) 



A 



n 

j=i 



a(A 



(JV-1) x(JV-2^ 



6(Af - x) - Af~ 2) ) 



(39) 



We would like to close this paper with the following remarks. First we note that it is 



possible to perform convenient shifts in the Bethe ansatz rapidities, {A^ } — ► {Af j } — |, in 
order to present the results (36-39) in a more symmetrical form. For instance, after these 
shifts, the nested Bethe ansatz equations can be compactly written as 

Lin 



A 



(a) 



'■'a . ir 

2 



\ ( a ) I Sa,w 



n n 



A 



(a) 



2 



\ (°) \0) I Ca,b 

6=1 k=l, k^i Aj — A fc H 2~ 



m a ; = 1, 



N — 1 



(40) 



where C a f, is the Cartan matrix of the Lie algebra and w — 1, N — 1. We note that for 
iV = 3 we recover the results by Abad and Rios ||. It is also straightforward to extend the 
above Bethe ansatz results to the superalgebra Sl(N\M). The N = case is the simplest one, 
since the only modification is the addition of minus signs in functions o(A) and 6(A). 

Next remark concerns the physical meaning of the ^-operators as scattering 5-matrices. 
It is known that C^i(X) might represent the scattering matrix of particles belonging to the 
fundamental representation of SU(N). The extra solution £^(A), however, should be seen as 
the forward scattering amplitude between a particle and an antiparticle. In fact, it is possible to 
show that the whole particle-antiparticle scattering (even backward amplitudes) can be closed in 
terms of the braid-monoid algebra. We leave a detailed analysis of this possibility, their Bethe 
ansatz properties as well as generalizations to include trigonometric solutions for a forthcoming 
paper [O 
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